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ABSTRACT 


Laser diffraction particle sizing was conducted at the nozzle entrance, inside the 
nozzle, and in the exhaust plume of a small two-dimensional rocket motor using two 
different metallized propellants to determine the effects of nozzle geometry on particle 
breakup. Six different nozzles were used, including three converging nozzles and three 
converging’diverging nozzles. 

An AP/GAP;AI propellant containing 4.69 percent aluminum showed no significant 
effects of nozzle geometry on exhaust plume particle size. Exhaust plume measurements 
showed consistent results under various conditions, indicating that the particle breakup 
had been completed at the nozzle throat. The particle size distribution was successfully 
tracked through the motor for one particular converging section, with particle breakup 
inside the nozzle occurring prior to that predicted by theoretical analysis. 

An AP. HTPB-ZrC propellant containing 1.0 percent zirconium carbide exhibited 
particle size variation in the exhaust plume, but sparse data and varying combustion 
chamber pressures precluded isolation of the causal factor. Additionally, combustion 
chamber pressure was shown to decrease measured particle size inside the motor for 
both propellants. 

Extensive particle sizing validation experiments using particles of known size sus- 
pended in distilled water were conducted throughout the course of the investigation to 


assist interpretation of laser diffraction particle sizing data. 
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I. INTRODUCTION 


With advantages such as simplicity, reliability, relatively low cost, and a well- 
developed technology base, solid propellant rocket motors are used in a wide variety of 
military and space vehicle propulsion systems [Ref. 1: p. 18]. Yet, solid propellants 
possess a limited specific impulse, prompting continued efforts to develop more energetic 
propellants. 

The introduction of metallic particles into solid propellants can provide significant 
increases in specific impulse. Figure 1 shows the theoretical increase in specific impulse 
from adding aluminum particles to a composite propellant [Ref. 1: p. 86]. Unfortunately, 
the resulting presence of metallic and metallic oxide particles in the combustion cham- 
ber, exhaust nozzle, and rocket motor plume creates additional concerns. 

In many cases, the particles perform the essential function of acoustic damping. 
Detrimental particulate effects include two-phase flow losses and increased visual and 
infra-red plume signatures. 

Acoustic damping minimizes the eflect of combustion oscillations inside the motor. 
These oscillations occur when two conditions are met. First, energv release mechanisms 
in the motor must amplifv disturbances. Second, the motor geometry and flow charac- 
terisucs must be capable of responding to these disturbances. The ensuing oscillations 
mav be either insignificant or so severe as to cause destruction of the rocket motor. 
Metallic (or non-metallic) particles suspended in the combustion chamber are capable 
of keeping motor operation stable in such conditions, provided that the frequency of 
oscillation is sensitive to the size and densitv of the particles present, and that a suflicient 
number of particles 1s available. [Ref. 1: p. 185] 

Two-phase flow losses are those losses caused by the presence of solid or liquid 
particles in combustion gases. Velocity lag is the largest contributor to two-phase flow 
losses: during expansion through the nozzle, particulates accelerate less rapidly than the 
Surrounding gas, resulting in a drag force on the expanding gas. Thermal lag also wastes 
available energy: particles maintain a higher temperature than the expanding gas, re- 
leasing a significant portion of their thermal energy after exiting the motor. Also, nozzle 
регіогтапсе 15 гедисед Бу (те inability of particles to expand. Effective means of esti- 
mating the above losses exist; the Air Force Astronautics Laboratory SPP (Solid 


Propellant Performance Prediction Code) yields loss estimates, but without a reliable 
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Figure 1. Effects of Aluminum on Theoretical Performance of AP/HTPB/AL 
Propellants 


means of deternuning particle size, a major source of uncertainty in. performance pre- 
diction remains [Ref. 2: p. 2]. 

Factors which affect the particle size distribution inside the motor are many, in- 
cluding: percent of metal in the propellant, metallic particle sizes m the unburned 
propellant, combustion chamber pressure, burning rate of the propellant, and residence 
time of combustion products in the motor. Exhaust particle size distributions may de- 
pend on the above factors, as wel] as the expansion process through the exhaust nozzle. 

In an attempt to determine the impact of the expansion process on particle size, ef- 
forts to size particulate matter in the combustion chamber, exhaust nozzle, and exhaust 
plume were undertaken using both converging and converging/diverging nozzles of 
varving geometry in a small “two-dimensional” rocket motor. Particle sizing was ac- 
complished at various chamber pressures using two sigmificantly different propellant 
compositions. Laser diffraction particle sizing was the primnary means of data 


acquisition. 


t 


The hostile rocket motor environment presents numerous challenges to effective 
particle sizing. Particle collection is one useful technique, but particle breakup during 
collection creates significant uncertainty. Additionally, physical access to the motor in- 
terior 1s required for sizing upstream of the exhaust nozzle. Laser diffraction particle 
sizing solves these problems, but not without compromise. The technique takes advan- 
tage of the strong influence of size on a particle's light scattering signature by illumi- 
nating a sample with laser light and measuring the diffraction pattern in the 
near-forward region. Particle sizing in motor exhausts by this technique is non-intrusive 
in that no particle breakup is induced, but interpretation of the data can be much more 
difficult than measuring the size of collected particles. Data can be taken inside a motor 
(virtually impossible using a particle collection technique) but the need for optical access 
presents further challenges. Windows must be placed in the motor and, much more 
significantly, these windows must be kept clean prior to and during data acquisition. 

A Malvern 2600HSD laser diffraction particle sizer was the primary means of par- 
ticle sizing used during this investigation. The Malvern uses a low power helium-neon 
(He-Ne) laser to illuminate a sample and measures the scattered light intensity as a 
function of scattering angle to enable calculation of a particle size distribution. Exten- 
sive efforts Were undertaken to ensure suitability of the Malvern system for particle 
measurements in the rocket motor environment. 

The primary objective of thesis research was to determine the effects of nozzle ge- 
ometrv on particle size distribution in the nozzle and exhaust of a small 2D motor. 
Secondary goals included evaluation of the Malvern 2600HSD system as an acceptable 
means of particle sizing in a rocket motor environment and determination of the effects 
of combustion chamber pressure and propellant composition on particle size 


distribution. 


H. BACKGROUND 


A. COMPOSITE PROPELLANT BEHAVIOR 
l. Aluminized Propellant 
Solid propellants are mixtures of fuel and oxidizer. The mixture may be on a 
molecular scale, as in homogeneous propellants, or on a much larger mechanical scale, 
as in heterogeneous composite propellants. Both types usually require additional ingre- 
dients to ensure acceptable physical properties. The propellant used most extensively in 
this investigation was heterogeneous GAP'AP/AL. [Ref. 1: p. 93] 


Table 1 shows the aluminized propellant composition. [Ref. 2: p. 24] 


Table 1. ALUMINIZED PROPELLANT INGREDIENTS 


mass 
Ingredient 
percent 


EGON 


ALUMINUM (C003) 
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Price [Ref. 3: p. 14-1] lists two advantages of including aluminum particles in 
heterogeneous propellants, as well as four disadvantages. Advantages are the potential 
for increased specific impulse and the suppression of acoustic oscillations. Disadvan- 
tages include erosion of motor components. poor combustion efficiency at low pressures, 
two-phase flow losses and an increased plume signature. Despite the disadvantages, 
aluminum is widely used in heterogeneous propellants. One further characteristic stated 
bv Price (and much of the literature) is that aluminum addition has only a minor influ- 
ence on propellant burning rate. 

GAP is both a binder and a fuel, containing large amounts of carbon and hv- 


drogen. Ammonium perchlorate (AP) 1s a “monopropellant” in that it contains both 


oxygen and fuel (but an excess of oxygen, so it 1s also an oxidizer.) Ammonium 
perchlorate is thus in itself a "fuel lean” propellant. Only a small amount of additional 
fuel 1s required to optimize the fuel'oxidizer ratio. The aluminized propellants used 
during testing contained more than 70 percent ammonium perchlorate by mass, making 
ammonium perchlorate combustion characteristics a major aspect of propellant burning 
behavior. 

Examination of burning ammonium perchlorate crystals by high speed photog- 
raphy and SEM analysis of quenched specimens has revealed the presence of a molten 
laver of AP on the burning surface [Ref. 4: p. 372]. This characteristic plavs an essential 
role in the agglomeration of aluminum particles on the propellant burning surface. 

Understanding the behavior of aluminum particles throughout the combustion 
process is essential im modeling two-phase flow losses and in estimating acoustic damp- 
ing capabilities. Many papers on the subject provide qualitative data, but no means of 
accurately predicting Al and Al,O, particle size in the combustion chamber, nozzle, and 
exhaust plume exists. Aluminum particle behavior 1s generally divided into three cate- 
gories: aluminum behavior on the burning surface, combustion during transit, and par- 
ticle breakup in the exhaust nozzle. 

a. Particle Agglomeration 

Caveny and Gany [Ref. 5: p. 13-1] describe a process by which aluminum 
particles remain in a reaction layer (the molten ammonium perchlorate) and form 
agglomerates, given three requisite conditions. First, particle diameters must be smaller 
than the reaction laver thickness. Second, surface tension forces associated with the re- 
action laver must be sufficient to overcome drag forces on the particles in the reaction 
laver. (Unburned particles must sustain motion through the molten AP to remain within 
the reaction layer. If the surface tension of the reaction layer is able to provide a force 
sufficient to overcome the resulting drag, particles will be carried along by the receding 
reaction layer.) Third, ignition of aluminum particles must be delayed. These three 
conditions can cause a significant number of aluminum particles to become trapped in- 
side the reaction layer. As the number of trapped particles increases, adhesion forces 
between the particles begin to form agglomerates. Oxide shells surrounding the particles 
can retard ignition, while cracks in the shells allow molten aluminum to flow between 
particles, strengthening the bond between them. Through this process large 
agglomerates may form, containing hundreds of the original aluminum particles. Large 
agglomerates are then ejected into the flow as their size exceeds the reaction laver 


thickness. While propellants behave uniformly in a macroscopic sense, the differences 


сл 


in AP crystal sizes and microscopic mixing anomalies combine to ensure a wide range 
of aluminum particle behavior on the propellant burning surface, from single particle 
ignition to the formation and ejection of large agglomerates. 
b. Combustion of Agglomerates 

Complete combustion of agglomerates in the combustion chamber produces 
sub-micron ALO, droplets which are not of great concern. Incomplete combustion of 
the agglomerates, however, can lead to relatively large particles which have significant 
impact on motor performance. “Burned-out” aluminum oxide shells can also contribute 
to flow losses. To ensure complete combustion of agglomerates in the combustion 
chamber, two conditions are necessary. First, the temperature must exceed 2300°K, the 
melting point of aluminum oxide. When this occurs, the oxide shell melts (if still intact) 
and retracts from the already molten aluminum agglomerate. Sufficient time in the 
motor (called residence time) is the second condition required for complete combustion. 


Price estimates the required residence time as: 
т = Ар" (1) 


In the above expression, t 1s the time required for complete combustion and 
D 1s the diameter of the agglomerate in microns. Typical values for the constants k and 


аге: 


п 1.75 
k = O.5x10^? sec [(um)" 


These constants predict a combustion time of 0.03 seconds for a 150 um 
agglomerate [Ref. 3: p. 13-4]. Aluminum particles, if burned completely, produce a mass 
of AZO, 1.89 times greater than that of the original aluminum particles. 

c. Breakup of Aluminum Particles 

Cavenv and Gany (Ref. 6: p. 1] used high speed photography through a 
small two dimensional rocket motor to record the breakup of aluminum and aluminum 
oxide agglomerates during nozzle transit. Experimental results verified the Weber 
number (We) as an important breakup predictor. 

The Weber number 15 a ratio of the inertial forces to the surface tension 
forces acting on a particle. The inertial forces result from a velocity differential between 


the particle and the surrounding gas. and induce particle breakup. The surface tension 


forces are a measure of the particle’s cohesion, or resistance to breakup. Calculation 
of the Weber number 1s shown below [Ref. 6: p. 2]: 


d, p, (u, — Ugg) 
И’е = ко з к (2) 


G 

where d,, — agglomerate diameter 
Pg = gas density 
u, = gas velocity 

H = agglomerate velocity 


o = agglomerate surface tension 


Cavenv and Gany found that particle breakup usually occurred when a 
Weber number of 28 was exceeded. While sub-micron particles can readily follow the 
flow, larger agglomerates typically breakup due to an inabilitv to accelerate rapidly 
enough to maintain a low Weber number. 

Even for large particles, the Weber number 1s typically low at the nozzle 
entrance. As the gas velocity increases, the particle lags. For nozzle geometries causing 
relatively gradual acceleration, Weber number 1s believed to accurately predict particle 
breakup. When acceleration of the gas 1s so rapid as to approximate exposure of the 
particle to an instantaneous gust, however, a high Weber number is insufficient to cause 


muumedrate breakup. A minimum breakup time, 7, , is needed [Ref. 6: p. 2]: 


Пота 





(3) 


While calculated breakup times are typically quite short, they are not in- 
significant When compared with nozzle residence times. The nozzle geometries used in 
the present investigation caused a rapid increase in Weber number near the throat. 
While it was believed that the high Weber number of larger particles ensured breakup, 
it was unclear what allowance, if any, for breakup time Was required. 

2. Zirconium Carbide Propellant 
Aluminum particles greatly increase the visible plume signature, making them 
inappropriate for use in some propellants. Even when a lower specific impulse is ac- 
ceptable, acoustic damping requirements dictate inclusion of particulates. Small 


amounts of zirconium carbide (ZrC) are often used in reduced smoke propellants. The 


second propellant used in particle sizing experiments was a ZrC composite. Major in- 


gredients of the propellant are given in Table 2 below [Ref. 7]: 


Table 2. MAJOR ZRC PROPELLANT INGREDIENTS 
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While similar in composition to the alununized propellant, the ZrC propellant 
had a significantly lower burning rate. The ZrC propellant contained a large amount 
of AP and was therefore characterized by the previously discussed "molten reaction 
laver. The longer burn times of the ZrC propellant probably allowed a greater per- 
centage of particles to remain in the reaction laver, but the propellant contained much 
less metal than the alununized propellant. The melting temperature of ZrC is signif- 
icantly higher than the adiabatic flame temperature for the propellant, indicating that 
the agglomeration process was more complex than for the alununized propellant. [Ref. 
7] 

ZrC dissolves in the presence of HCl and HINO, [Ref. 8: p. 24]. While 
ammonium perchlorate (NH,ClO,) has potential to form both. chemical equilibrium 
calculations [Ref. 9] indicated that HNO, was not an expected product of the ZrC 
propellant combustion. Netzer [Ref. 7] states that the most likely agglomeration se- 
quence for ZrC commences with the liberation of oxvgen from the ammonium 


perchlorate during combustion. The oxygen then reacts with ZrC as follows: 


ZrC + O, => Zr + CO, 
or (3) 
РИС чє О, -2 745 -- 220 


Oxygen then further reacts with the exposed zirconium: 


Zr +O, > ZrO, (4) 


The complete reaction, shown below, has an adiabatic flame temperature of 
approximately 4400°C [Ref. 7]: 


BC oO Eo" CO, (5) 


The adiabatic flame temperature for the propellant calculated by Micropep [Ref. 
9] was 3018°K. For comparison, Netzer lists the melting and boiling temperatures 


shown in Table 3. 


Table 3. ZRC ALLOY MELTING AND BOILING TEMPERATURES 
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The melting temperatures of Zr and ZrO, are both sufficiently low to ensure that 
they exist in a molten state in the reaction laver. It was not believed that ZrC was di- 
rectly involved in the agglomeration process due to its high melting temperature. Parti- 
cle sizing experiments verified significant agglomeration in the ZrC propellant, and it 


was believed that the agglomerates Were primarily composed of ZrO, . 


В. LIGHT SCATTERING 

The manner in which particles scatter light depends upon characteristics of the par- 
ticle, the illununating light, and the suspending medium. Mie’s theory predicts light in- 
tensity as a function of scattering angle, size, internal structure, refractive index (relative 
to the medium), and incident wave polarization and wavelength [Ref. 10: p. 1]. Figure 
2 [Ref. 11: p. 2] shows attenuation of incident light on a spherical particle bv scattering 
and absorption. 

When the wavelength of the illuminating light and the diameter of the scattering 
particle are approximately equal, simplification of Mie theory can lead to unsatisfactorv 
results. With a particle diameter significantly different than light wavelength. however, 
simplification of Mie theory provides excellent results. Most noteworthy is that forward 


scattered light at small angles 1s accurately described by Fraunhofer diffraction if the 
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Figure 2. Attenuation of Light by a Spherical Particle [Ref. 11 : p. 2] 


particle diameter is significantly larger than illuminating light wavelength. In this case, 
the refractive index of the particle has no effect on light scattering. The pattern of light 
scattered by the particle becomes a strong function of particle diameter and the wave- 
length of the incident light, while essentially independent of other factors. The scattering 
signature of a particle changes rapidly with particle size in this regime; the measurement 
of forward scattered light at small angles can therefore provide an excellent means of 
particle sizing s ЕС p | 

For polvdispersed particles illuminated by. coherent monochromate hght, the re- 
lationship between scattered light intensity, scattering angle, and particle diameter 15 as 


follows IRe TE p0): 





со 2 
1, DIE | 2 (20) | 
— = = СОЕ 0 
: ës - (o) (6) 


where /, = intensity at angle 0 


1) = intensity of original unscattered light 
0 = forward scattering angle 
a = rad 


= particle diameter 
д = illuminating light wavelength 
Jı = Bessel Function of order 1 


N(a) — the particle size distribution 


Equation 6 requires independent scattering, an acceptable approximation when 
particles are separated by a distance of at least three radi [Ref. 13: p.5]. Additionally, 
single scattering is assumed. When particles are densely packed іп а sample, multiple 
scattering effects become significant. The extent to which multiple scattering effects can 


be removed from collected data 1s discussed in the Experimental Apparatus section. 


C. SIZE DISTRIBUTION PARAMETERS 

Comparison of particle size distributions 1s simphfied if standardized, descriptive, 
easily obtained parameters are used. For a completely identified distribution (1.е. the 
exact size of each particle is known) utility becomes the single most important factor in 
the selection process. For example, Hermsen |Кеї. 14: p. 2] maintains that D,, . the 
mass-weighted average diameter, 1s the most useful size distribution parameter for pre- 
dicting two-phase flow losses. D,, 1s calculated as shown below: 
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D is the particle diameter and n is the total number of particles. As an example, 
D,; for a distribution containing 20 particles with a diameter of 4.0 microns, 15 particles 
with a diameter of 8.0 microns, and 5 particles at 12.0 microns would be calculated as 


follows: 


= 29.67да (5) 


11 


This value ts considerably larger than the average particle diameter, which is 6.5 microns. 
D,; 18 a mass Weighted parameter and will always be larger than the average diameter for 
any distribution except one in which all particles are of identical size. 

Calculation of such parameters is a simple matter provided that the distribution is 
known. This is often not the case. Many previous sizing efforts did not obtain the size 
distribution explicitly. For example, Cramer et al. [Ref. 15: p. 13] discuss two techniques 
whereby a particle size distribution parameter 1s obtained without prior determination 
of the distribution. In both cases, the parameter obtained was D,, , the Sauter mean 
diameter (SMD). Cashdollar et al. [Ref. 16: p. 1764] describe a similar means of ob- 


taining D, : 


С! 
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where f = transmittance 


О = average extinction cociiicient 
C, — mass concentration of particles 
{ = path length 


p = particle density 


Р; — Sauter Mean Diameter 


Cashdollar states that Q can be closely approximated knowing only D,,. Thus, if 
С, and p are known, transmittance can be measured over a path length to obtain a value 
for D4. While РД, тау provide a better means of predicting two phase flow losses, Ds, 
is more prevalent in the literature dwe to its experimental accessibility. D,, is calculated 
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D, for the size distribution of the previous example ts calculated as follows: 


20(-4.0)? + 15(8.0)? + 5(12.0) 
О. = o X c DE 8.00um (11) 
СОР Шр 


D. 16 less than D4, but still considerablv larger than the average diameter. Phys- 
ically, D,, 1s six times greater than the ratio of the total volume to the total surface area 


for a set of particles. 


ПІ. EXPERIMENTAL APPARATUS 


A. THE MALVERN PARTICLE SIZER 
1. General Description 
a. Optical Components 
The Malvern 2600HSD laser diffraction particle sizer makes exclusive use 
of the Fraunhofer diffraction simplifications to Mie scattering [IRef. 17: p. 17]. The sys- 
tem consists of an optical transmitter, an optical receiver, a synchronization unit, and a 
desk top computer with associated software. Figure 3 [Ref. 18: p. 1-11] depicts the op- 


tical sub systems. 
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Figure 3. The Optical Measurement Unit [Ref. 18: p. I-11] 


The optical transmitter and optical receiver combine to form the optical 


measurement unit, specifications for which are given below [Ref. 18: p. 1-15]: 


e Laser transmitter. 2.0 mW He-Ne laser (633 nin wavelength) with 9.0 mm beam 
expander 


e Receiver. Fourier transform lens mount, detector on sliding rail carriage, and de- 
tector electronics 


e Rangelens. Interchangeable Fourier transform lenses, three available focal lengths: 
63 mm, 100 mm, 300 mm 


e Detector. 3] element solid state detector array, optimum concentric semi-circular 
annuli design 


e Detector electronics. 31 amplifier parallel sample’ hold construction, A/D conver- 
sion and on board digital storage, 10 ws sampling ume for all rings, 25 ms data іп- 
put time, sampling from either internal or external trigger source 

b. Principles of Operation 
Fraunhofer diffraction dominates in the forward scattered lobe throughout 
the applicable size limits of the Malvern particle sizer. Within the limits of this as- 
sumption, the variation of light intensity with scattering angle is suflicient to determine 


a particle size distribution. Figure 4 demonstrates this concept. 
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Figure 4. Light Scattering by Nlonodispersed Particles 


Illuminated particles scatter light in all directions. Large particles, however, 
tend to scatter a large percentage of the forward scattered light at small angles. Con- 
versely, smaller particles scatter more light at relatively high angles. In addition to il- 
lustrating this concept, Figure 4 shows how the Malvern svstem exploits. this 
phenomenon. The Fourier transform lens forms tlie far field diffraction pattern of the 


scattered light at its focal plane [Ref 18: p. 1-7]. Simply put, it converts scattering angle 


to radial displacement. All light reaching the Fourier transform lens will be focused on 
the detection array with a radial position determined solely by its incident angle to the 
lens. Light scattered at the same angle from different scattering centers will impinge 
upon the detection array with identical radial displacement from the beam axis. An ar- 
ray of 31 concentric semi-circular detection rings 1s centered about the beam axis on the 
lens focal plane. A measure of the light energy striking each detector ring thus gives 
scattering intensity as a function of 31 sets of scattering angles. 

Figure 5 shows Malvern presentations of scattering data for two 
monomodal distributions of polystyrene spheres suspended іп distilled water. Note that 
for the smaller particles, the signature shifts towards larger angles (the outer rings). 

The manner in which collected data is converted to a particle size distrib- 
uuion 1s proprietary and not available from Malvern. A Malvern application note states 
that the distribution of scattered light 1s compared with calculated distributions for var- 
lous size distributions until an optimum fit 1s obtained [Ref 19: p. 1]. The effect of this 
process 15 to invert /(@) versus (@) data to produce 31 size bands from the measured 
scattering aneles [Nem 1k: 242 |: 

c. The Fourier Transform Lens 

Three lenses were available for use with the 2600ITSD model, each with 

different particle sizing abilities. Table 4 lists the lenses by focal length and indicates the 


size ranees Or edach. 


Table 4. LENS CHARACTERISTICS 


Focal Sieg Cut-off Approximate Nlaximum 
m W Range Distance = Angle 
| 63mm | mm 1.2 to 118 microns 55 mm | tsa a 











300 mm 


Proper lens selection is critical to ensure that dynamic range is optimized 
without violating cut-off distance. The shorter the lens focal length, the smaller the 
particles that can be measured. This sensitivitv to smaller particles is due to the ability 
to detect scattered light at higher angles. Unfortunately, the measurement of a wider 
range of scattering angles necessitates that the scattering centers be closer to the lens to 


ensure that all measurable scattered lieht strikes the lens. Sample placement outside the 
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Figure 5. Scattering Signature for Two Different Monomodal Distributions 


cut-olf distance results in bias towards larger particles because the light scattered at high 
angles (bv small particles) misses the lens. The maximum scattering angle measured by 
cach lens was not provided by the manufacturer. Values shown in table 4 are estimates, 


obtained trigonometrically as shown in Figure 6. 


lens diameter 
beam diameter 

= lens cut-off length 
maximum scattering angle 


278 


0 ~tan! 


max 


Figure 6. Maximum Scattering Angle Calculation 


To ensure collection of scattered light at the maximum collection angle for 
all particles in the control volume, it 1s sufficient that particles at the intersection of the 
cut-off distance and the upper edge of the beam are visible to the range lens at the 
maximum scattering angle. The exposed diameter of each lens was used in conjunction 
with manufacturer provided cut-off distances to calculate an approXinate maxmum 
scattering angle measurable bv each lens. 

d. Malvern Operation 

Several basic steps arc common to all particle sizing experiments. These 
include selection of the appropriate lens, specification of user-defined parameters, system 
alignment, background measurement, sample measurement, and size distribution calcu- 
lations. After selecting the appropriate range lens using the previously discussed criteria, 
user-defined parameters are set. These parameters include experiment tvpe (whether 
measuring particles in air or water) active beam length (the length of beam inside the 
sample), and distribution type. Previous laser diffraction particle sizers required user 
specification of a distribution type (i.e. log normal or Rosin Ramler.) The Malvern 
26001150 has a model independent option which enables sizing of multimodal distrib- 
utions. 

With the system in the intended configuration, an X-Y translator positions 


the detection array such that the incident beam is centered on (Пе атау. Гсало а 


15 





small variations in the beam path due to in-line optics. Once the system 1s aligned, a 
background reading is taken. The measured light intensity readings for each detector 
ring are stored for later subtraction from the sample measurement data. A sample 
measurement can then be taken under the desired conditions. 

Sample measurement can be triggered internally (by a user at the svstem 
keyboard) or externally (bv any desired condition, such as a motor chamber pressure 
threshold). During sample measurement the Malvern averages light intensity values on 
each ring over a user-defined number of sweeps. Array sweep time is approximately 10 
us, but processing time increases this value to about 7 ms. The number of sweeps taken 
in any experiment is a trade-off between statistical significance and data acquisition 
duration. 

Subtraction of the background levels from each averaged ring intensity 
value produces intensity versus scattering angle data. Under ideal conditions, the re- 
sulting profile is a unique function of particle size distribution alone. System software 
then converts light intensity versus scattering angle data to a volume distribution and 
provides several size distribution parameters. Figures 7 and 8 show results of key steps 
in the process. Figure 7 shows a background measurement and a sample measurement. 
Figure 8 shows the calculated size distribution data and a volume distribution plot. 

2. System Limitations 
Typically, the Malvern 2600HSD laser diffraction particle sizer is used for 
droplet sizing in spray patterns, or the measurement of particles in liquid. The 2D rocket 
motor environment creates some difficulties for the Malvern system, requiring careful 
analvsis of all data to remove bias and prevent over-reliance on Malvern generated re- 
sults. The most significant difficulties encountered are discussed as follows. 
а. High Obscuration 

Obscuration represents the fraction of light which is scattered or absorbed 

by particles in the beam path. Transmittance is the fraction of light which is not scat- 


tered or absorbed. The two are related as follows: 
Obscuration = 1 — Transmittance (12) 


Malvern software divides sample measurement center diode intensity by the 


background value to obtain transmittance: 


T Ee Sample Measurement Center Diode Intensitv (13) 
ЕНСЕ —————— s 
Background Measurement Center Diode Intensity 
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Malvern Background Measurement 
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Malvern Sample Measurement 
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Figure 7. Data Acquisition 


The Malvern measurement is reported in terms of obscuration, which is an 
excellent indicator of the amount of multiple scattering occurring in a sample. If the 
obscuration is low, for example 0.05 or five percent, then 95 percent of the transmitted 
light is unscattered. Of the five percent which is scattered, some is reflected or absorbed 
by particles and some is diffracted. Much less than five percent of the original light is 


actually diffracted in the forward lobe, striking the range lens. As 95 percent of the or- 
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The Calculated Size Distribution 
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Figure 8. Calculated Distribution Presentation 


iginal light passes through the sample unscattered, it is reasonable to expect that the 
light which is scattered once has a similar or higher probability of passing through the 
remainder of the sample without being scattered a second time. This single scattering 
condition allows accurate Malvern particle sizing, and exists in practicality up to an 
obscuration of 50 percent. Ideal Malvern measurement obscurations are between five 


percent and SO percent--below five percent there are an insufficient number of scattering 


2] 


centers In the sample volume, and above 50 percent multiple scattering effects begin to 
become significant. 

The multiple scattering which appears at the higher obscurations causes 
calculated particle size distributions to be shifted towards smaller particles. As a signif- 
icant amount of scattered light becomes rescattered, the light intensity signature shifts 
from the inner rings to the outer rings of the detection array. Malvern software makes 
no corrections for high obscuration. Changing the obscuration value for a data set (by 
inputting a false obscuration value but leaving all other values the same) does not 
change the particle size distribution calculations. This is an acceptable condition only 
if the user does not violate the 50 percent obscuration limit. 

Exhaust gas particle sizing experiments Were not adversely affected by mul- 
tiple scattering, as obscurations during such runs were routinely within the ideal range 
or just over 50 percent. Data taken through the motor, however, often had obscuration 
values near 99 percent. Such data were more qualitative than quantitative, but still al- 
lowed important comparisons. Empirical corrections for high obscurations are discussed 
in the Validating Experiments section. 

b. Sensitivity to Combustion Gas 

The sensitivity of Malvern hardware to exhaust gas produced by the 2D 
motor is not a shortcoming of the equipment (the stainless steel motor was not 
Impervious to the corrosive effects of HCl in the exhaust products) but nevertheless it 
Was a condition requiring attention. The entire optical system was encased in an alu- 
minum box with windows placed over optical access holes to prevent exhaust gas from 
entering. Two disadvantages of this configuration were the decrease in available cut-off 
range and obscuration of the optical path by additional windows. For data acquisition 
through the motor, the beam passed not onlv through densely packed particles, but also 
through four windows (two on the motor and two on the protective enclosure.) While 
background measurements compensated for the windows, less power was available for 
sample measurement. More significantly, additional optical elements increased the 
likelihood of multiple reflections. Alignment of the system was more difficult, and the 
possibility of secondary reflections degrading data acquisition accuracy was increased. 

с. Sensitivity to Combustion Light 

The detector array was not capable of discriminating between He-Ne light 
and light near the He-Ne wavelength. The system was sensitive to sunlight (the test cell 
door was completely open for motor firings) and to combustion light produced by the 


motor. The effects of sunlight were minor due to the Malvern orientation in the test cell, 


to 
n 


and were easily eliminated by proper background measurements. Motor light was a 
more difficult problem because it was inconsistent and inseparable from the effects of 
light scattering by particles. A wide-pass filter was therefore incorporated to block mo- 
tor light. A narrow pass filter was not used as transmittance of the filter would have 
been a strong function of incidence angle, grossly biasing results. Use of the wide-pass 
filter was a trade-off between passing a small amount of other than He-Ne hght and bi- 
ased light attenuation as a function of scattering angle. 

Figure 9 shows a comparison between exhaust runs with and without the 
filter under similar conditions. Figure 10 shows the effects of sunhght on the system 
with and without the filter. Data acquisition without the filter showed effects similar to 
those obtained by exposing the system to sunlight: an increase in outer ring intensity 
readings. The similarity in these results lead to the conclusion that combustion light 
grossly affected unfiltered data acquisition. Filter use during the sizing of particles in 
liquid verified that the filter did not change the scattering signature for non-luminous 
particles. Hence, the filter was added to the inline optics, bringing the number of optical 
elements between laser and range lens to five. 

d. Large Minimum Measurable Particle Size 

The wavelength of the He-Ne laser is 0.6328 um. Recalling that the 
Fraunhofer diffraction simplification to Mie scattering theory 1s valid only for particles 
with diameters significantly larger than illuminating light wavelength, successful sizing 
of particles less than one micron would not be anticipated. Table 4 confirms this ex- 
pectation, showing the smallest sizeable particle to be 1.2 microns for the 63 mm lens 
and 1.9 microns for the 100 mm lens. Malvern logic provides for a sub-class of particles 
and estimates the amount of particulate matter below the minimum measurable size that 
IS present in the sample. Such calculations are inexact, however, and can lead to un- 
certainty in particle sizing results. Typical particle sizing attempts in the 2D motor 
showed a trimodal distribution of particles, with the smallest mode in the sub-micron 
range. The abilitv of the particle sizing unit to handle such small particles is further 
discussed in the Validating Experiments section. 

e. Short Cut-off Distances 

Wablew4 sows tlie culoll distances for the range lenses. The 63 mm lens 
cut-off distance is 55 mm, or roughly 2.2 inches. Motor geometry and the required 
Malvern protective enclosure made it impossible to present the entire sample volume 
within 2.2 inches of the range lens. Attempts at data acquisition with a slight violation 


of the 63 mm lens cut-off distance were frustrated by ambiguous results. Data taken in 
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Figure 9. Combustion Light Filtering 


this configuration showed a rapid decrease in light intensity with scattering angle at high 
scattering angles. 

Figure 1] shows exhaust run results under similar conditions for both the 
63 mm lens and the 100 mm lens. Data taken using the 63 mm lens shows a decrease 
in intensity with scattering angle at the outer rings. While it 1s possible that the pub- 


lished cut-off distances provided a “safety-factor” and that the intensity versus scattering 
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Background Measurement Without Filter 
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Figure 10. Sunlight Filtering 


angle distribution using the 63 mm lens was correct, it is more probable that excessive 
distance from the range lens caused the light scattered at high angles to miss the lens. 
The figure also shows the effects of beam steering on the inner rings. This phenomenon 
I C liSsed separately. 

The lack of ambiguity in the 100 mm lens data was deemed to outweigh the 


benefit of the smaller sizing capability of the 63 mm lens. All test data were therefore 
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Exhaust Data with the 100 mm Lens 
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Figure 11. Exhaust Data Using Different Lenses 


collected using the 100 mm lens, with the penalty of a larger minimum measurable par- 
ticle size. 
f. Wide Beam Width 
The 9 mm He-Ne beam used by the Malvern system may have combined 


with motor geometry to provide a slight bias during runs through the motor. Addi- 
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tionally, the wide beam width made the beam placement decision a difficult one during 
exhaust runs. 

Video taped motor firings revealed the presence of normal shocks in the 
exhaust plume for all nozzle geometries. Typical distance between the motor exit plane 
and the first shock was 0.4” or just over 10 mm. 

Placement of the 9 mm beam as close as possible to the aft end of the motor 
would have caused some scattered light to be blocked or reflected by the motor. Pro- 
viding sufficient separation between motor and beam to alleviate this condition meant 
that a portion of the sample volume would be affected by the shock wave. Beam 
placement in exhaust runs which were intended to measure particle size as close to the 
nozzle exit as possible was therefore a compromise between the two adverse conditions. 
Data at the exit plane was taken with the beam positioned between 0.06" and 0.41” be- 
hind the motor. 

Inside the motor, effective optical access (clean windows) required that the 
free boundary between the pressure chamber and the window cavity be as small as pos- 
sible. The motor was designed with a 9 mm diameter free boundary between the pres- 
sure chamber and each window cavity. Light scattered at the edges of the beam, if 
scattered away from the beam axis, was blocked at the entrance to the far window cavity 
after passing through the combustion chamber. With the 8.4° maximum collection angle 
of the 100 mm lens and a 0.25” beam length, the bias towards larger particles bv this 
effect was believed acceptable. 19.4 percent of the sample volume was restricted in some 
way by this condition. Half of this amount (the bottom half) was restricted only in di- 
rections that would not strike the detection array. Also, the estimated maximum meas- 
urable scattering angle for the lens is an absolute maximum, and is probably high. Some 
hght scattered at lower angles was also blocked, reducing the bias towards larger 
particles. 

g. Beam Steering 

All scattered light reaching the range lens is used in particle size calcu- 
lations. While most of this light is scattered by diffraction, some scattering 15 the result 
of refraction due to large density gradients between the exhaust plume and the sur- 
rounding atmosphere. All exhaust runs produced intense scattering on the inner rings 
which rapidly decreased with scattering angle. This phenomenon has been identified by 
Malvern, and is due to the beam being steered off-axis by densitv gradients at the rela- 


tively narrow free expansion boundaries. Figures 9 and 1] show the impact of beam 


steering on exhaust data. The sharp decrease in intensity with scattering angle near the 
beam axis is not a characteristic of any particle distribution scattering signature. 

Effects of beam steering were reduced by the “kill data” command, which 
allowed the system to calculate a particle size distribution without considering light in- 
tensity values for a given number of either the inner-most or outer-most rings. The kill 
data command killed ring data in pairs, up to a maximum of ten rings on either end of 
the array. For exhaust runs, the kill data command was used to ignore data on the inner 
rings affected by beam steering. Reliance on outer ring scattering data was not expected 
to affect data interpretation due to the small size (and large associated scattering angles) 
of exhaust particles. 

Unfortunately, not all effects of beam steering were removed by the kill data 
command. Diflracted light was also subject to beam steering, and it 1s believed that the 
effects were similar to those of high obscuration values--a shifting of the scattering in- 
tensity towards higher angles. As a result, beam steering was expected to shift the cal- 
culated particle size distribution towards smaller particles. An additional effect of beam 
steering is an artificially high obscuration, produced by knocking some of the otherwise 
unscattered light off axis. While no quantitative data was available, it was believed that 
correction of the particle size distribution for the artificially high obscuration caused by 
beam steering indirectly corrected for the beam steering bias. 

3. Laser Diffraction Validation Experiments 
а.  Мопотоаа! Distribution Sizing Capabilities 

Preliminary Malvern calibration was accomplished using monomodal dis- 
tributions of polystyrene microspheres in distilled water. Previous Malvern users had 
raised doubts concerning the applicability of laser diffraction techniques to transparent 
particles, but svstem performance in sizing the particles was excellent. 

Polystyrene microspheres suspended in liquid were obtained from the Duke 
Scientific Corporation. The particle distributions were of two tvpes. The first type was 
monomodal, with particle diameters varying several microns from the mode. These 
particles were documented by Duke with a histogram indicating the relative number of 
particles in each of approximately twenty size bands. Figure 12 shows the distribution 
histogram for 10.4 micron particles. The second distribution type contained particles 
of much more uniform size. These distributions were documented with a mean diameter 
and a standard deviation. The particles used in monomodal validation experiments were 


of both tvpes and ranged in size from 20.0 microns to 0.109 microns. 


XX 

XX 

XXX 

XXXXX 

KXXXXXX 

KXXXXXXXX 
KXXXXXXXXXXX 
KXXXXXXXXXXXXXX 
XXXXKKXXXKXXXXXXXX 
XXXXXXXXXXKXXXXXXXXX 
XXXXXXXXXXKXXXXXXXXX 
XXKKXXXXXXXXKXXXXXXX 
XXXXXXXXKXKXKXXXXXXX 
XXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXXX 
XXXXEXKKXXXXXX 
XXXXKXXXXXXX 
XXXXXKXXX 

XKXKXXXX 

XXXXXX 

XXXX 

XXX 

XX 

X 

X 





anm 


ор 
el 


to 


Ф 


i 


15 


button И 


str 


cle D 


art 


on P 


i 


C] 


10.4 NI 


е 12. 


jgu 


1 


\ 


v SE 
stribution of 


ee hy 


lon 


al distribut 


SCV eT 


r 


О 


f 


Won та | 


leu correct 


f 


Cri 
s Duke document 


aise Cie У 


COST 


Duke h 


- 


ih wide di 


а 


r 
с 


Figure 13 confirm 


photography. 


[ 5.00 


ОШО 


stribut 


о лош ог М! 


5 


Figure 14 


particles about a 10.4 nucron mean. 


meron particles. 


Эс 


NAST 


"47121217, 





29 


SEM Analysis of 10.4 Micron Nlicrospheres 
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Figure 14. SENI Analysis of 35.00 Micron Müicrospheres 


Sizing was simplified through use of the “PS1” &пттеб ратпсе пе аш 
The cell. designed for use with the Malvern system. has an electrically driven magnetic 
stirrer and is easily mounted on the optical measurement unit. For monomodal distrib- 
uuon experiments, the cell was filled with distilled water and placed in the beam path. 
After svstem alignment, suspended particles were added to the cell one drop at a time 
until reaching an optimum obscurauon. The sample was then measured and a distrib- 
ution calculated. 

Figure 15 shows a plot of Malvern calculated D,, versus actual D,, for the 
polvstvrene spheres using the 100 mm lens. The results are excellent for particles larger 
than the published 1.9 micron minimum measurable size. For smaller particles, Malvern 
“sub-class” logic appeared to consistently assign a D,, value of 0.9 microns to the 


distribution. 
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Figure 15. Malvern Monomodal Sizing Perforinance 


b. Trimodal Distribution Sizing Capabilities 

With the prevalence of bimodal and trimodal distributions in past particle 
sizing experiments, it was decided to evaluate Malvern capability to handle such dis- 
tributions using particles in liquid. Before receiving the PSI particle stirring unit, ex- 
periments were run with miniature sample cells constructed of microscope slides and 1/4” 
square strips of Plexiglas. Sinular experiments were later conducted using the PSI stirred 
cell. Scanning electron microscopy was also used to further validate results and to ex- 
amine Malvern sizing performance using aluminum particles. 

(1) Miniature Cell Trimodal Experiments. Trimodal distribution exper- 
Iments using the miniature cells were conducted by creating a known trimodal mix and 
using 1t to evaluate Malvern sizing ability. For cach experiment, three monomodal dis- 
tributions of polystyrene particles in distilled water were sized in 1.5 ml sampling cells. 
S Iis ofeach cell were then poured into a larger 4.5 ml cell, creating a trimodal 


nux for which a size distribution was determined by the Malvern. Results were com- 
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pared to an anticipated value of D,, based upon Malvern obscuration, manufacturer 
provided size distributions for the polystyrene particles, and a calculated scattering effi- 
ciency for the particles. The governing theoretical expression is discussed Dy Gomi IRER 
20: p. 3555] and is identical to Equation 9: 


E 2 —3 0 С,! 
I bes exp( 2 р,» ) (14) 


where ///, = transmittance 
= Sauter Mean Diameter 
= optical path length 


average extinction coefficient 


Oo. $ 
|| 


— volume concentration of particles 


For each monomodal distribution of polystyrene parucles, all values 
in the above expression except C, were known or calculable. The necessary valuessies 
transmittance, optical path length, and D,, in Equation 14 were readily obtained or 
measured; calculation of the extinction coefficient was much more complex. However, 
once these values were determined, Equation 14 vielded C,. The manufacturer provided 
size distribution, the cell volume, and the volume concentration of particles were then 
sufficient to describe the distribution. The theoretical trimodal size distribution was then 
easilv calculated bv combining the three monomodal distributions. 

The optical path length, /, was equal to the internal dimension of the 
miniature cell and was easily measured. Sauter Mean Diameter was calculated nem 
manufactured supplied data for each set of particles (and verified bv SEM for some dis- 
tributions.) Transnittance was determined from obscuration measurements made bv the 
particle sizer. While using Malvern results to evaluate Malvern performance was unde- 
sirable, the obscuration calculation is a simple one involving only two measurements of 
beam intensity, easily accomplished by the system. Malvern obscuration measurements 
Were Validated by measuring the obscuration through filters of known optical density (at 
632.8 nm.) Figure 16 shows Malvern transmittance values for several of the filters, and 
includes a curve representing the theoretical transmittance as a function of the optical 
density. 

Calculation of Q, the remaining unknown in Equation 14, was more 


involved. Q is the average extinction coefficient for a distribution of particles. The ex- 
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Figure 16. Validation of Malvern Obscuration Values 


tinction coefficient of a single particle is simply Q, and is the cross-sectional arca of the 
light scattered or absorbed by a particle in terms of its flat plate area. A particle which 
absorbs all incident light with no effect on light which passes near it would have a value 
Гог О of 1.0. Diffraction of nearby light allows for extinction coefficient values greater 
than one. In fact, large particles typically have cxtinction coellicients of 2.0. 
Cashdollar Ref, 16; p. 1764] desernibes the means by which Q 1s cal- 
culated. Mic scattering theory (not simplified diffraction theory) and numerical inte- 
gration yield Q for a particular particle size given the relative refractive index. of the 
particle and the wavelength of illuminating light. The refractive index of the polystyrene 
microspheres is 1.587 at 540 nm. Although the Ile-Ne wavelength is 632.8 mm, use of 
the refractive index at 540 nm was believed to provide acceptable accuracy [Ref. 21]. 
There is no imaginary component for the index of refraction of the particles because the 


particles are transparent. 


Miescat [Ref. 22], a Fortran program written by Cashdollar and re- 
vised by Weldon, allows calculation of O for an assumed log normal distribution given 
the refractive index of the particles relative to the medium. The program was further 
modified to allow calculation of Q as a function of particle size and Q for a user defined 
particle size distribution. Figure 17 shows the scattering efficieney of a ро 5 ШШЕ 
sphere as a function of its diameter for Ile-Ne lrght as calculated by the revised Miuescat 


program. 


= 
раї 
Еа 
— 
є» 
м 
Gs 
f=, 
E 
- 
O 
7 
C 
— 
T 
© 
A 
— 
E 
У 
Е 


7.5 10125 POSEE MEER 
DIAMETER IN MICRONS 





Figure 17. Polystyrene Microsphere Scattering Efficiency 


Cashdollar gives the expression for numerical calculation of Q over 


a number of size bands shown below [Ref 16: р. 1764): 


3 O(d) N(d)d^ Ad 
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where Q(d) = particle extinction coefficient 
N(d) 2 number of particles in size band 
d = particle diameter 
Ad = size band width 


For each monomodal distribution of particles, the manufacturer's size 
distribution and the index of refraction were used to calculate D,, and Q as previously 


described. Table 5 lists the calculated results for each monomodal set of particles. 


Table 5.5 MONOMODAL POLYSTYRENE SIZE DISTRIBUTIONS 
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The size distribution volume listed in the table is the total volume of 
all particles listed in the manufacturer's size distribution. For the smaller particles, 
which were documented only by a standard deviation, a size distribution was generated 
from the standard deviation. For each monomodal distribution, volume concentration 
values from Equation 14 were multiplied by the cell volume to obtain total particle vol- 
ume. The number of particles in each size range of the manufacturer's distribution was 
multiplied by the ratio of the total particle volume to the size distribution volume to 
obtain the correct number of particles in each size band for each monomodal distrib- 
ution. Calculated distributions were then combined to obtain a theoretical trimodal 
distribution. from which D; was calculated. The following tables show the results of 


several experiments and allow comparison of Malvern and theoretical data. 
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Table 6. TRIMODAL MIX NUMBER ONE 


Distribution 
Malvern Volume 


Di; (um) Transmittance Concentration 


08295 925103 


ҮШҮ 02510 
19.6 1557.10: 


Malvern Mixture D,, - 6.5шп 
Theoretical Mixture D, = 11.1 um 





Table 7. TRIMODAL MIX NUMBER TWO 


Distribution 
Malvern Volume 


Dz, (um) Transmittance Concentration 


0.624 0.1166 ].198x10-* 


Malvern Mixture D,, 2 2.4 um 
Theoretical Mixture D,, — 1.9 um 





Table 8S. TRIMODAL MIX NUMBER THREE 


Distribution 
Malvern Volume 


Di; (um) Transmittance Concentration 


Malvern Mixture D,, 6.2 шт 
Theoretical Mixture D,, 2.6 шт 





Table 9. TRIMODAL MIN NUMBER FOUR 


Distribution Malvern Volume 


Ds, (um) Transmittance Concentration 


— 029 | 630910 


5 0.6791 1.034х10-* 
0.6936 1762x107 


Malvern Mixture D,, 
Theoretical Mixture D,, 





The trimodal mixture shown in Table 6 was created without violating 
either the dynamic range of the lens or the recommended obscuration limits. Even 
though the mixture obscuration was just above five percent, the Malvern calculated D, 
was lower than expected. Multimodal distributions such as the one depicted in Table 6 
tended to reduce calculated values of D,, , and did so in a repeatable fashion. Figure 18 
shows Malvern size distributions for each of the monomodal distributions of Table 6 as 
well as the Malvern trimodal distribution. The shaded peaks in the plot are the indi- 


vidual monomodal distributions. 





Figure 18. Comparison of Monomodal and Trimodal Results 
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Improper mixing of the particles could have affected calculated values 
of Dj, , but would not have moved the mode peaks. The systematic underestimation of 
mode peaks in trimodal distributions was consistently repeated, and lead to underesti- 
mation of D, for multimodal distributions. Nevertheless, the correct identification of 
the trimodal distribution type was most significant. 

In sinular experiments, Gulder [Ref. 23: p. 2] showed that for bimodal 
distributions, increased obscuration values tended to shift the calculated volume per- 
centage of particles from the larger mode to the smaller mode. Thus, Malvern calcu- 
lations Were not expected to give precise Values for the individual mode peak diameters 
or the exact volume percentage associated with each mode, but the relative relationship 
between mode peaks was believed accurate. 

Tables 7, $, and 9 show the effects of reducing the smallest mode size 
below the minimum measurable particle size. This models actual conditions in the 2D 
motor. Table 7 depicts results for a mixture containing a 0.624 um mode. While sig- 
nificantlv below the minimum measurable particle size for the lens, the sub-micron par- 
ticles clearlv affected the measured D,,. Malvern and theoreucal values for Doer m 
and 1.9 microns, respectively. These values are close enough to be the result of exper 
imental error, but the fact that the Malvern calculated D,, was higher than the theore- 
tical value is significant. Had the Malvern been able to accurately measure the 
sub-micron particies, the Malvern generated D,, would have been lower than the the- 
oretical value (due to multimodal effects). Ignoring the sub-micron particles, the the- 
oretical D,, for the remaining bimodal distribution was calculated at 15.5 microns. Thus, 
While the Malvern did not accurately measure the smaller particles, Malvern results 
much more closely approximated the true distribution than an approximate distribution 
containing onlv the larger particles. 

Tables $ and 9 show interesting results obtained using particles of 
approximately 0.5 microns in diameter and smaller. Such small particles were unde- 
tected bv the Malvern in multimodal distributions. The effect was more pronounced 
with the 0.109 micron particles of the distribution in Table 9. The Malvern D,, of 6.2 
microns was what would have been expected for a multimodal mix if the smaller particles 
had not been present. The theoretical D,, for the distribution considering only the larger 
particles was 7.4 microns. The Malvern value was slightly lower than this, as would have 


been expected duc to the multimodal effects on calculated D,,. 
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Figure 19 shows a Malvern generated volume distribution plot for the 
distribution of Table 9, which appears to be a bimodal distribution of the larger particles. 
The Malvern did not detect the smaller particles. For such small particles, much of the 
scattered light misses the range lens, making a slightly higher obscuration the only 


discernable contribution to the light scattering signature. 





Figure 19. . Trimodal Distribution with 0.110 Micron Particles 


Figure 20 shows Malvern ability to detect the smallest mode in a 
multimodal distribution as a function of its size. Each value of Malvern SMD is plotted 
(as an “x”) against the smallest mode diameter in the trimodal mix. The vertical lines 
on the plot represent the corresponding theoretical D,, calculations. The bottom point 
of each line represents the actual theoretical D,, of the trimodal mix. The top point 
corresponds to what would have been the D,, if the smallest mode were not present. (It 
1s the theoretical bimodal SMD for the two larger modes.) While the volume percentage 
of the smallest mode varied from distribution to distribution, smallest mode size effects 
are clearly shown: if the smallest mode is comprised of 0.5 micron particles or smaller, 


it does not affect multimodal measurements. 
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Figure 20. Malvern SMD for Distributions with Small Modes 


(2) Stirred Cell Trimodal Lxperbnents. lhe increased mixture control 
of the PSI stirred cell was believed significant enough to warrant further trimodal sizing 
efforts. The procedures used were similar to those for the miniature cells. Slight mod- 
ification Was required, however, as only one stirred cell was available i 
monomodal distributions suspended in equal volumes of distilled water was sized in the 
stirred cell. The three monomodal distributions were then nuxed to provide three iden- 
tical trimodal distributions. All three trimodal mixtures were sized to ensure consistent 
nuxing. The theoretical results were calculated as before. A scanning electron micro- 
scope was then used to validate the theoretical results for the trimodal mix. Table 10 
shows the results for a trimodal mixture consisting of 20.0, 9.6, and 4.8 micron particles 


and reconfirms the previously noted node dependent decreases an 


Table 10. STIRRED CELL TRIMODAL RESULIS 












И Malvern Theoretical 
Particle 
Distribution D32 Ds, 














[le sir еШ оаа е were dehydrated and sized using a 
scanning electron microscope. Comparison of SEM results with theoretical and Malvern 
calculations further validated the theoretical approach and confirmed that Malvern sizes 
Were consistently small for trimodal distributions. Figure 21 shows one of three photo- 
graphs taken of the trimodal mix. 

The three particle size modes are clearly visible in Figure 21. The 
number of particles in each mode was counted for three similar photographs and an ac- 
tual D,, was calculated. For calculation purposes, all particles in each distribution were 
estimated to have a diameter equal to the mean particle diameter. Table 11 shows the 


particle count. 


Table 11. SEM TRIMODAL RESULTS 


Mean Particle 
Diameter 
(microns) 


Number of 
Particles 





A total of 814 particles were counted. It would have been desirable 
to count a greater number of particles for more statistical significance, but the resulting 
Dj, of 13.53 microns agreed exceptionallv well with the theoretical 13.77 microns. Thus, 


the stirred cell SEM experiment not only provided further evidence of the relationship 
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Figure 21.  Trimodal Distribution SEM Photograph 


between calculated D,, and the number of modes, but also validated the transmittance 
equation used for obtaining theoretical data. 
c. Dispersed Distribution Sizing Ability 

A sample of aluminum particles idenucal to those used in casting the 
aluminized propellant was sized using the sturred cell. Nfalvern Sizing results were then 
compared with SEM photographs to evaluate Malvern performance. ligure 22 shows 
the distribution calculated bv the Malvern. 

The Malvern calculated a dispersed distribution extending from 6 microns 
to the lens limit of 118 microns with a D,, of 29.2 microns. Figure 23 SHom coom 
photograph of a portion of the same sample that was sized by the Malvern. Three such 
photographs were taken. Each particle in the three photographs was measured and a 


D;, for the distribution was calculated. 
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Figure 22. Malvern Size Distribution of Propellant Aluminum 
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Figure 23. SEM Photograph of Propellant Aluminum 
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The photographs showed only 414 particles, with diameters between two 
and 60 microns. Most particles appeared to be in the five to 20 micron range. 50 of the 
particles were larger than 20 microns, producing a D,, of 23.74 microns. This compared 
reasonablv well with Malvern calculations. 

d. Obscuration Corrections 

Obscuration values obtained through the motor were well in excess of those 
recommended bv the svstem manufacturer. An empirical correction for the multiple 
scattering effects of high obscuration developed by Gulder [Ref. 23: p. 1] was evaluated 
using the stirred cell and various distributions of spherical particles. 

Gulder states that the correction is applicable for trimodal as well as 
monomodal distributions. Gulder obtained the following expressions by curve fitting 


numerous Malvern sizing results: 
F, + F 
C, = 13524 tD) 


13.1226 5.74740 
Di; ) D ) 
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D os 
F, = -—041184(—) + ( 


100 uu 


Р, = 2.2289: - 2.60779” 


In the above, $ 1s the percent obscuration and D,, 1s the value calculated 
бу the Malvern. Adjusted D,, 1s obtained by dividing the Malvern D,, by C, , the cor- 
rection factor. The correction 1s applicable for obscurations between 50 and 98 percent 
and for Malvern measured D,, values above ten microns. 

Figure 24 shows the relative magnitude of the correction. The ratio of the 
actual D,, to the reported D, is plotted for reported Sauter mean diameters of ten and 
30 microns. Throughout the range of application, the Gulder correction indicates actual 
values for D,, between one and two times ereater than Ще пером a mi 

Figure 25 depicts results of a Gulder correction validation experiment. 32.0 
micron aluminum oxide particles supplied by the K.C. Abrasive Company were sized by 
the Malvern at high obscurations using the stirred cell. The figure shows reported Dj, 
and corrected D,, as functions of obscuration. Results were somewhat disappointing. 

Other such experiments using polystyrene particles produced much better 
results. It was believed that the poor results obtained using aluminum oxide particles 


were caused by dynamic conditions during sample measurement. The aluminum oxide 
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parucles were exceptionally difficult to keep suspended in distilled water, even with 
benefit of the magnetic stirrer. Ilence, obscurations measured bv the Malvern (as well 
as other data) were averages of significantly different individual readings. Unfortunately, 
similar dynamic conditions were believed to exist inside the motor. 

]he Gulder correction provided estimates of actual D,, that were better than 
“order of magnitude”, but still inexact. The most serious restriction to its application 
was the required minimum reported D,, of 10 microns. Attempts at extending the range 


of application proved unsuccessful. 
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Figure 24. Obscuration Effects on Reported SMD 
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Figure 25. Gulder Corrected SMD for Aluminum Oxide Particles 


e. Beam Steering Compensation 

Effects of the kill data command were investigated by comparing Malvern 
results with and without the command using the same raw data. Only the smaller dis- 
tributions were examined because the beam steering phenomenon (and thus the need to 
kill data) was observed only during exhaust runs. 

Figure 26 shows the use of "Kill 5,07, which caused svstem software to ig- 
nore the first ten rings. In this case. use of the command increased the calculated D,, 
from 10.2 microns to 10.4 microns and shghtly changed the mode shape. 

Smaller particle distributions were unaffected by the КШ data command 
Additionally. no remarkable effects were noted in multimodal distributions which had 


al} mode sizes below 2U microns. 
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Figure 26. Use of the Kill Data Command 


f. Validating Experiment Conclusions 

Lessons learned during experiments with particles in liquid had direct ap- 
plieatron to 2D motor particle sizing experiments. Treh obscuration and the presence 
Of multimodal distributions were shown to cause a decrease in measured values of D,,, 
With the obscuration effects being more sigutficant. Iligh measured values of D,, were 
obtained when the smallest mode was less than approximately 0.5 microns. The kil] data 
command was shown to have little impact on the sizing of small particles, making the 
major eflect of beim steering (high readings on the inner rings) relatively 
inconsequential. 

cece мое expected to Бс шо! аПссісі by 
obscuration. Reasonable compensation for the obscuration was believed possible, but 
ШИ спе 101105 of Gulders correction. Further. accurate indication of all size 
 атгеутпап 025 пустой was anticipated, but the specific Malvern mode sizes were 


treated as approximations. 


B. THE TWO-DIMENSIONAL MOTOR 

The two-dimensional rocket motor was designed by Pruitt [Ref 17: p. 28] and most 
recently used m conjunction with laser diffraction particle sizing by Hovland [Ref. 2]. 
The motor was used to burn two 1°4” by 3/4" propellant slabs of various lengths and 


allowed data collection internally as well as in the exhaust plume. A nitrogen purge 
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system maintained optical access (clean windows) during runs through the motor. No 
purge was required for exhaust data collection. Motor ignition was accomplished using 
small BKNO, igniters installed at the head end of the motor. Six sets of nozzles were 
used interchangeably, each consisting of two 1/4" flat plates held in place bv stainless 
steel posts. On all runs, the potential for motor over-pressurization was reduced by in- 
stallation of a disposable 1000 psi burst disk. 

Figure 27 shows the motor interior with nozzle plates installed. The 1/4” deep motor 
cavity was 11.6 inches long from head end to exhaust plane. Motor height was 2.5 


inches at the head end and 2.3 inches near the nozzle. 





Figure 27. The Motor Interior 


The two holes visible in Figure 27 along the longitudinal axis of the combustion 
chamber allowed installation of a splitter plate to ensure propellant ignition bv directing 
igniter discharge towards the propellant surfaces. Use of the splitter plate was found to 
be necessary only for runs which did not use nitrogen purge. 

Also visible in the figure are two 4.3 inch propellant slabs. During actual motor 
loading. propellant surfaces in contact with motor walls were coated with commercially 


available RTV high temperature gasket maker to provide a burning inhibitor and to 
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maintain propellant position in the motor. Surfaces not touching motor walls were left 
uncoated. A regressive burn resulted from leaving the aft slab ends uncoated, but coat- 
ing the slab ends would have increased the likelihood of introducing RTV into the 
combustion products. 

1. Igniters 

Motor igniters were constructed from 7/16-20UNF bolts, insulated copper wire, 
nickel cadmium wire, and BKNO,. A 3/16" diameter cavity was cleared inside the shaft 
of each igniter bolt and two 0.031” holes were drilled through the bolt head into the shaft 
cavity. An insulated copper wire was inserted through each 0.031” hole from the bolt 
head into the shaft cavity and a short (approximately 1/8”) length of nickel cadmium 
wire was soldered between the copper wire ends. The wires were then positioned such 
that the nickel cadmium igniter filament was centered in the shaft cavity and secured by 
applying epoxy glue to the bolt head exterior. Igniter construction was completed by 
filling with BKNO, and glueing a paper cap over the igniter cavity. Motor firing was 
accomplished by bolting the igniter to the head end of the motor and placing a 12V DC 
source across the copper wires. Heating of the nickel cadmium resulted in combustion 
of the BKNO, and subsequent propellant ignition. 

2. Nitrogen Purge 

The nitrogen purge system was used during internal runs to prevent window 
fouling. Figure 28 shows the window assembly components. The bottom pieces are, 
from left to right. a tapered window cavity insert, a cylindrical sintered bronze nitrogen 
diffuser, a window receptacle, and a cover plate. For data collection in the motor ex- 
haust, the window cavity was blocked with the plug shown at the top of the figure. (No 
purge was used for exhaust measurements.) 

Figure 29 shows both motor halves. The motor half containing the nozzle 
plates and the combustion chamber cavity is plugged, while an unplugged window cavity 
can be seen in the other motor section. 

Each window assembly received high pressure nitrogen through a 0.031” sonic 
choke. The nitrogen flow was diffused through the cylindrical sintered bronze insert and 
exited the window cavity through the tapered window insert. Nitrogen supply pressure 
was set at 1200 psi (greater than twice expected chamber pressure) to ensure constant 


flow through the chokes during motor operation. 
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Figure 28. The Window Assembly 
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Піешге 29. The Disassembled Nlotor 
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Hovland calculated the nitrogen mass flow rate for each window using the 
choked flow formula with high pressure nitrogen at room temperature and a circular 
sonic choke throat diameter of 0.031" [Ref. 2 : p. 20]: 


қ ШЖ; E (y -1)/(y—1) ы 
mom Pg po num (17) 


where P, = nitrogen pressure 
А, - sonic choke area 
у = 1.4 
ge = 32.174 ft Ib, /Ib, sec 
7, = nitrogen total temperature 
К = 55.0 ft lbglb,, R° 


The calculated nitrogen flow was 0.02 ]b,/ sec per window. Similar calculations 
for a combustion chamber pressure of 500 psi using the aluminized propellant indicated 
a combustion gas flow rate of 0.13 lb,/ sec. Thus, the ratio of purge gas to combustion 
gas for internal runs was approximately 0.3. This high purge flow rate probably pro- 
vided signifrcant quenching of the particles. This was not considered a problem since 
no measurements were made further downstream when purge was used. 

3. Nozzle design 

A single converging ‘diverging nozzle designed to provide cross-sectional area 
changes similar to that of a 3D 30° converging nozzle was designed by Pruitt [Ref. 17: 
p. 30] and used by Hovland. Lessons learned from this nozzle provided the basis for 
subsequent nozzle designs. Specific problems encountered using the original nozzle in- 
cluded a large uncertainty in expected chamber pressure and difficulty interpreting data 
taken through the motor due to the sharp nozzle curvature in the sample volume. 

The inability to accurately predict chamber pressure led to burst disk ruptures 
on several occasions, invalidating data. Such occurrences were believed due in part to 
the small nozzle throat area. The relationship between throat area and chamber pres- 
sure 1S: 

Ақару“ 1/(1—n) 


(18) 


Eig 


"| 


where P, — chamber pressure 

A, = burning surface area 

Au = nozzle throat area 

Ыр = propellant density 

g. = gravitational constant 
с = characteristic exhaust velocity 
a — burning rate constant 


n — burning rate pressure exponent 


Equation 18 shows that as nozzle throat area is increased, errors in measuring 
propellant slab length (and thus A,) become less significant. Increasing the nozzle throat 
area was also believed to reduce the possibility of momentary nozzle obstruction. 
Equation 1$ was used to calculate the maximum throat area which would permit the 
desired range of chamber pressures. Maximum propellant slab length (and therefore 
A,) was limited by motor geometry. The largest possible throat area calculated to allow 
chamber pressures of $00 psi (considering the limitations on burning surface area) was 
0.04 square inches (0.16" x 0.257.) 

The original nozzle is shown in Figure 30. The drastic mitial converging angle 
was required to maintain the same cross-sectional area versus length characteristics as 
found in a converging 30° three dimensional nozzle. 

The original nozzle was deemed unsuitable for studying the eflecr piki 
geometrv on particle size. Sharp turning at the nozzle entrance modeled cross-sectional 
area characteristics of a three dimensional nozzle, but not the flow patterns. It was not 
known to what extent particle breakup in the nozzle was induced by the sharp turns. 
Furthermore, cross-sectional area varied greatly over the 9 mm beam width, and the 
nozzle entrance cut through the sample volume. While this still allowed investigation of 
the effects of chamber pressure and propellant composition on particle size, it was of 
little value in studying the effects of nozzle geometry on particle size. 

Nozzle design was conducted with the goal of accurate particle sizing upstream 
of the nozzle, inside the nozzle, and in the motor exhaust. Constant slope converging 
sections were deemed to be best suited for initial investigations for several reasons. First, 
anv measured particle breakup would less hkely be caused by sharp turning of stream 


lines or by agglomerates striking the nozzle plates. Second, for data taken inside the 





Figure 30. The Original 2D Nozzle 


nozzle, conditions across the sample volume were much more uniform. Linally, the 
straight converging sections simplified nozzle fabrication, allowing more time for data 
acquisition. 

The three converging nozzles designed and built for use in the 2D motor are 
shown in Figure 31. Each had a throat area of 0.04 square inches and a 0.050 inch flat 
at the throat to nunimize effects of nozzle erosion. 

A preliminary analysis using combustion chamber stagnation conditions and 
one-dunensional idealized expansion (1.c., perfect gas, adiabatic, and single phase) was 
undertaken in an attempt to predict particle breakup in the nozzles. The short and long 
converging sections were compared. The first step in the simplified model was the cal- 
culation of gas velocity as a function of position inside the nozzles. 

Zucker (Ref. 24: p. 120] derives expressions describing steady, adiabatic, once- 
dimensional, isentropic flow of a perfect gas. Gas velocity is the product of Mach 


number and acoustic velocity. Mach number was found from the following: 
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Figure 31. 
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The Converging Nozzles 





[26-0 (15) 


A(x) = Area as a function of position 
A” = nozzle throat area 


Mi = о нее 


= rauo ol specific l мк 


x = distance from nozzle entrance 


For each x, A/A* was determined from nozzle geometry. The ratio of specific 
heats, y , was calculated by Micropep for the aluminized propellant and assumed con- 
stant throughout the expansion process. An iterative scheme produced Mach number 
as a function of position. 


Calculation of the acoustic velocity, a, was also facilitated by Micropep output. 
a = Jg, RT (20) 


R is the gas constant and T 1s the temperature of the gas. The gas constant and 
combustion chamber temperature were provided by Micropep. The combustion cham- 
ber temperature was used as the stagnation temperature, 7, , to determine temperature 
throughout the nozzle: 

1, 


r= ——ə  . .. (21 
ро юю — 1 


Gas velocity inside the nozzle was simply the product of Mach number and 
acoustic velocity. Particle velocity (as a function of size) was then calculated using gas 
velocity and a predictor-corrector tvpe scheme in conjunction with an approximate ex- 
pression for particle acceleration in terms of the gas/particle velocity differential. 


Caveny and Gany [Ref. 6: p. 3] approximate the drag force on a particle as: 
F = 0.5Cpp,(V, — V,)°(a/4)d; (22) 


where p, — gas density — P/RT 


Й о = gas velocity 

Є 1 > Lea 
F, = particle velocity 
q, = particle diameter 


dhe drag coellicient, C, , 1s dependent upon the Reynolds number, R, . With 


particle diameter as the characteristic length, the Reynolds number is defined: 


ST eer А 
К coe (23) 


where и = gas viscosity 
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Cavenv and Gany found the drag coefficient as follows: 


Ср = 27184 (ог В, 380 
= 0.271R,*" Tor 10 5 EO (24) 
= 2.0 for R,, > 10° 


Similar approaches have been taken bv other authors. Rudinger [Ref. 25: p. 


122] gives an identical expression in terms of particle acceleration: 


Es 
б ж as n s ay (25) 
D, 4а, Ор 75 р 

A computer program was written taking the above approach to determine both 
gas velocity and particle velocity as functions of position in the nozzle. Weber number 
could then be found as a function of position and particle size for each nozzle. (See 
Equation 2.) Figure 32 shows the Weber number as a function of position for three 
particle sizes in the short and long converging nozzles. 

The curve representing the 1.5 micron particle is not visible in the figure: it was 
found to follow the flow, with a negligible Weber number at the nozzle exit. The 20 and 
40 micron particles were unable to do so, with high exit Weber numbers resulting. The 
figure shows a similarity between the nozzles: significant particle (and gas) acceleration 
Is restricted to the throat region. Figure 33 further illustrates this point by showing the 
one-dimensional gas velocitv inside each nozzle. 

Similar calculations for particle diameters ranging from two to 50 microns 
produced the plot of throat Weber number versus particle diameter seen in Figure 34. 
Сахепу апа Gany observed breakup at a Weber number of 28, although breakup was 
expected between Weber numbers of 12 and 20. The horizontal line in Figure 34 corre- 
sponds to a Weber number of 28. It was believed that shear forces were sufficient to 
cause breakup of any particles exceeding this value. Particles with diameters larger than 
12.7 microns in the short nozzle would exceed this critical Weber number, as would 15.7 
micron particles in the long nozzle. 

Anticipated breakup characteristics based on throat Weber number were there- 
fore similar for all three nozzles. Such was the case for other designs considered as well. 
One significant difference among the three nozzles was particle residence time. The 
straight converging sections Were expected to allow examination of associated effects bv 


comparison of results for nozzles with significantly different residence umes. In addi- 
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Figure 32. Weber Number Through the Converging Nozzles 


tion, meaningful particle size measurement inside the nozzle would not have been pos- 
sible for other nozzle designs considered. For these reasons, the original plans to use 
Straight converging sections were carried out. 

Alter mal sizing experiments with the converging nozzles, two 
converging/diverging nozzles were designed to further investigate the effects of nozzle 
gcometr¥y on particle breakup. The medium length converging nozzle was chosen as the 
basis for comparison. Two converging/diverging nozzles with the same converging sec- 
tion and a 15° half angle expansion were designed and built, one with a converging sec- 
tion producing ideal expansion to atmospheric pressure for a combustion chamber 
pressure of 150 psi and the other with a design pressure of 500 psi. Both nozzles are 


enum ligure 35. 
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Figure 33. One-Dimensional Gas Velocity Through the Converging Nozzles 
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Figure 34. Throat Weber Number as a Function of Particle Size 





Figure 35. The Converging/Diverging Nozzles 
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IV. DATA ACQUISITION 


Unusual hardware used in this investigation was limited to the previously discussed 
2D motor and Malvern 2600HSD laser diffraction particle sizer. Less remarkable but 
nevertheless essential supporting equipment which formed the data acquisition network 
is discussed in this section, as well as normal data acquisition procedures. 

Figure 36 shows the data acquisition signal flow and all associated hardware. Each 
component in the figure is discussed in the following paragraphs in the order in which 


it was used during the data acquisition process. 


A. PREFIRING PROCEDURES 

Prefiring procedures commenced with calibration of the Teledvne pressure 
transducer. Pressure transducer output voltage at ambient conditions was amplified by 
a Pacific model 8655 amplifier and measured by two independent components. First, the 
amplifier voltage was measured and recorded using a Doric Integrating Microvoltmeter. 
Amplifier output was also recorded at the Hewlett-Packard 3054A Automatic Data 
Acquisition Control Svstem console, which continuously monitored amplifier output. 
After recording ambient readings, a 500 psi pressure was applied to the transducer using 
a dead weight tester. Amplifier output was again recorded and a linear slope intercept 
calibration was calculated for both monitoring devices. 

The second step in the prefiring procedure was the insertion of updated calibration 
data into pressure monitoring software. MALI, an HP basic computer program, was 
run on the Hewlett-Packard 9836 computer with the purpose of triggering Malvern 
measurements during desirable combustion chamber conditions. Ihe updated 
voltage pressure conversion was inserted into the code to ensure accurate measurement 
of chamber pressure. During a run, MALI monitored pressure transducer output 
through the Data Acquisition’Control System and triggered data acquisition at a user 
specified delav time after motor pressure reached a desired threshold value. 

Voltmeter calibration data was used to update Labtech Notebook, a commercially 
available data acquisition package used on one of the IBM PC AT computers. Labtech 
Notebook recorded amplifier output over the entire run, facilitating pressure-time trace 
construction. 

Proper configuration of Malvern accessories was the next step in prefiring proce- 


dures. The signal generator was set to produce a 100 KIIz triangular waveform and the 
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Figure 36. Data Acquisition Signal Flow 


Malvern Spray Synchronizer was set at 0.01 ms. This allowed a user specified number 
of externally triggered sample measurement sweeps with a minimum delay between 
sweeps. 

After clearing the motor firing area, the loaded motor was mounted on a test stand 
between the Malvern transmitter and receiving optics. For exhaust runs, the Malvern 
optical unit was positioned for desired beam location in the motor exhaust. For runs 
through the motor, the Malvern was positioned to ensure an unobstructed beam path 
through the motor windows with no reflected light reaching the diode array, and the high 
pressure nitrogen system was connected to both window cavities. In either case, system 
alignment was then accomplished using the detection array X-Y translator to compen- 
sate for inline optics. 

A background measurement completed Malvern preparation. Malvern software, 
which was installed in an IBM PC AT, measured the light intensity on each of the de- 
tector rings and saved the results for subtraction froin sample measurement values. For 
runs through the motor, the purge system was activated during background 
measurement. 

The prefiring sequence was completed by bolting an igniter to the motor head. Prior 
to doing so, the surrounding area was verified clear for safety considerations. Zero 
voltage across the ignition leads was confirmed using a digital voltmeter. The ignition 


leads were then connected to the igniter. 


B. THE FIRING SEQUENCE 

The motor firing sequence began with activation of a safety alarm. MALI was then 
run on the HP Automatic Data Acquisition Control Svstem. The program monitored 
transducer pressure in a continuous loop, awaiting the specified Malvern trigger condi- 
tions. Labtech Notebook data acquisition was then initiated. The program recorded 
pressure as a function of time for a 20 second period. The Malvern “Trigger External, 
Measure Sample" commands were executed, causing the Malvern to monitor Sprav 
Synchronizer output for data acquisition cues. Fire control panel power was activated 
and, for internal runs, nitrogen purge was initiated. The remote fire button was then 
pushed, placing a 12 V DC voltage across the igniter leads. 

Data acquisition during the run progressed automatically, based upon the specified 
trigger pressure and delav time in the MALI code. Successful propellant ignition caused 
combustion chamber pressure to exceed the specified threshold pressure (if specified 


correctly.) MALI then sent trigger signals to the Sprav Svnchronizer Unit. In addition, 
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a simultaneously generated reference spike was monitored by Labtech Notebook to re- 
cord the exact time of data acquisition. Malvern data was collected over the specified 
number of sweeps (typically 50), completing the process. 

The HP 2631B printer produced a real time printout of chamber pressure. More 
useful pressure-time data was provided by Labtech Notebook, which also indicated data 
acquisition time. Figure (37) shows a typical pressure time trace for alumnum 
propellant plotted from Labtech Notebook data. Initiation and termination of Malvern 


measurements are denoted with a “+”. 
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Figure 37. Typical Pressure Time Trace 


V. RESULTS AND DISCUSSION 


During the course of this investigation, 54 motor firings were conducted. This in- 
cluded 38 runs which vielded usable sizing data for either the ZrC or aluminized 
propellant. Five of the remaining measurement experiments produced useful sizing in- 
formation, but test conditions were not comparable to those of the primary data set. 
The results of two of these runs are included in the Experimental Apparatus section, one 
showing the effects of 63 mm lens use and the other showing the detrimental impact of 
unfiltered combustion light on raw data. The other three of these experiments were 
conducted with propellants other than the ZrC or aluminized propellants. Five motor 
firings were conducted without the Malvern to allow recording of exhaust plume phe- 
nomena with a video camera. Data acquisition for the remaining six motor firings was 


unsuccessful. 


A. ALUMINIZED PROPELLANT RESULTS. 

The study of aluminized propellant behavior was conducted in four parts: Malvern 
and SEM sizing of aluminum particles identical to those used in propellant casting, ex- 
amination of combustion bomb data obtained by Arnold [Ref. 26], internal 2D motor 
sizing, and exhaust plume measurements. This approach provided a means of “follow- 
ing” the particle size distribution through the combustion process, showing sizes of: 

e Particles in the unburned propellant 

¢ Quenched particles, collected immediately after leaving the burning surface 
ө Particles upstream of the exhaust nozzle 

e Particles inside the nozzle 


ө Particles in the exhaust plume 


l. Unburned Particles 
A discussion of the preburned particle size appears on page 42 in conjunction 
with Malvern dispersed distribution sizing capability. SEM and Malvern sizing 
produced Sauter mean diameters of 23.7 and 29.2 microns, respectively. An SEM count 
of 414 of the particles showed none larger than 60 microns, while the Malvern showed 
the upper limit of the distribution to be just over 100 microns. The previous discussion 
referred to a Malvern volume distribution plot of the particles. A Malvern number dis- 


tribution plot (obtained by mathematical transformation of the volume distribution plot) 


is shown in Figure 38. Figure 39 shows a similar plot constructed using SEM results. 
Comparison of the number distribution plots shows agreement between the two sizing 


techniques, increasing confidence in the data. 
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Figure 38. Malvern Number Distribution of Propellant Aluminum 


2. Combustion Bomb Data 
Arnold [Ref. 26] burned strands of the aluminized propellant in a combustion 
bomb to determine the size distribution of particles leaving the propellant burning sur- 
face. The strands were mounted on an inverted post inside a chamber which was pres- 
surized with nitrogen. Ignition was accomplished with a nickel cadmium filament. 
Burning particles were simultaneously collected and quenched in distilled water imme- 
diately beneath the burning surface. Table 12 shows SEM sizing results for the 


aluminized propellant at two pressures. 


Table 12. COMBUSTION BOMB RESULTS 






The 500 psi D,, was obtained from 814 counted particles. The 250 psi value was 
produced by only 59 particles, far fewer than deemed necessarv for statistical signif- 


icance. The low D,, values could have been caused by incomplete quenching or particle 
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Figure 39. Number Distribution of Propellant Aluminum Using SEM 


breakup during collection. Such was not believed to be the case due to reasonable 
agreement with Malvern data collected upstream of the exhaust nozzle. Arnold's data 
indicated that there was no significant agglomeration on the surface of the aluminized 
propellant strands. 
3. Internal 2D Motor Results 
A breakdown of test conditions and results for internal particle size measure- 
ments is shown in Table 13. Ten successful particle size measurements were conducted 


through the motor using the aluminized propellant. 
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Table 13. INTERNAL ALUNIINIZED PROPELLANT RESULTS 
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Run number eight in Table 13 pointed towards the high dependence of 
agglomeration on chamber pressure during tail-off. The Malvern was accidentally set 
Гог 200 sweeps rather than the normal 50, producing a data acquisition period of roughly 
1.3 seconds. Consequently. a significant amount of data were collected during 
propellant burn-out. The high measured D, was believed due to increased 
agglomeration at the lower pressure. The results of this particular size measurement 
were not considered in nozzle comparisons. 

Data taken using the short converging nozzle (convergent section length: 0.55 
inches) Were measured upstream of the nozzle. It was believed that this data could ac- 
curately be substituted for entrance conditions to other nozzles due to similarities in 
nozzle throat area (and therefore entrance Mach numbers.) Differences 1n motor resi- 
dence times for the different nozzles were originallv accounted for by the use of spacers 
inside the motor. Spacer use was seen to have no effect on size data; however dandas 
discontinued when experiments with the converging diverging nozzles began. Data col- 
lected upstream of the short converging nozzle were therefore used as nozzle entrance 
data for the medium converging nozzle. 

The converging diverging nozzle designed for expansion to atmospheric pressure 
from a chamber pressure of 150 psi had a converging section identical to that of the 


medium converging nozzle (convergent section length: 1.05 inches.) For the reasons 


С» 
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just discussed, it was believed that beam position was the only difference between 
measurements taken through the motor with either nozzle. Thus, 150 ря 
converging ‘diverging nozzle data were used as converging nozzle data with a sample 
measurement location corresponding to beam position in the 150 ря 
converging/diverging nozzle. 

Figure 40 shows the variation of particle size with position for internal data 
collected using the short converging nozzle, the medium converging nozzle, and the 150 
psi converging/diverging nozzle in terms of the corresponding measurement position for 


the medium converging nozzle. In the figure, data are presented as follows: 


р,, ° F. ‚ Obs. 










Ë а 
artíiala aise (am). 


D,, 1s the Sauter mean diameter, P, 1s the average combustion chamber pressure 
during data acquisition, and Obs. is the measured sample obscuration. Pressure and 
obscuration are included due to their impact on D,,: high pressure causes a real de- 
crease in D,, by restricting agglomeration and high obscuration causes a decrease in 
measured D,, due to multiple scattering effects. A miniaturized volume distribution plot 
is also included with each data set. 

The volume distribution plot gives the volume percentage of particles in each 
size band. Figure 8 on page 21 shows the range of each band for the 100 mm lens. The 
ordinate scale is automatically set by the Malvern, with maximum values of either 20 
or 100 percent. Unfortunately, both scales were used to display 2D motor data. 

Five of the six measurements shown in Figure 40 were conducted at average 
chamber pressures in the 350-400 psi range. One measurement was conducted at 627 
psi, with a correspondingly low D,,. The high pressure was most likely caused by RTV 
inhibitor failure or excessively thick (one or two hundredths of an inch) propellant slabs. 

Three measurements were conducted upstream of the nozzle (1.e., using the 
short converging nozzle.) One had an obscuration of 0.97. The other two obscurations 
were 0.995, too high to produce reliable results. As expected, the measurements with 


high obscurations produced notably smaller values of D,,. The remaining three meas- 
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Figure 40. Aluminized Propellant Internal Data; Medium Converging Nozzle 
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urements depicted in Figure 40 were conducted under similar conditions, but at different 
positions. The breakup process inside the nozzle was clearly shown. both in the D,, de- 
crease and in the shift of the volume distribution plot towards smaller particles. The 


breakup sequence is repeated for clarity: 


ЕШ 255 .97 











1% іт 
Particle sine (am), 





The measured values of D,, underestimate the true particle size. The 
obscurations were within limits of the empirical Gulder correction, but measured particle 
sizes Were not. No attempt to determine an exact correction for obscuration was made. 
Experience with the Gulder correction indicated that the measured values of D,, were 
roughly half the actual values. 

It was not apparent from the volume distribution plot why D,, was only 3.7 
microns for the upstream case. (It appeared from the distribution that it should be at 
least 20.) The volume distribution plot did not clearly show the volume of particles be- 
low 1.9 microns. For this measurement, the sub size volume (the volume of particles 
below 1.9 microns) was 24.6 percent. Figure 41 shows a number distribution plot for the 
Same data, which indicated that while the sub size particles comprised only 24.6 percent 


of the volume, they accounted for the vast majority of the number of particles. 
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Figure 41. Number Distribution of Upstream Particle Size Measurement 


The volume percentage of particles in the sub size for measurements proceeding 
towards the nozzle throat were 33.8 and 53.2, respectively. The volume distribution 
plots showed that this increase in the sub size volume was the result of agglomerate 
breakup: the volume distribution plot shifted towards smaller particles as the throat was 
approached. 

Breakup inside the nozzle was not predicted by theoretically obtained Weber 
numbers for either of the beam positions in the nozzle. In fact, methods discussed in the 
Experimental Apparatus section produced a Weber number of 2.0 for a 100 micron 
particle at the beam position closest to the throat. Such a low Weber number would 
not be expected to cause breakup (or even distortion.) 

Figure 42 shows results of the lone measurement inside the long converging 
nozzle (convergent section length: 2.0 inches.) It was believed that the low D,, was the 
result of high chamber pressure rather than breakup in the nozzle. 

The results in Figure 43 show measured particle size distributions using the ori- 
ginal nozzle, which simulated a conical converging section. While not directly compa- 
rable to the data collected in the straight converging section nozzles, the measurements 


verified the impact of pressure on Dy. 
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Figure 42. — Aluminized Propellant Internal Data; Long Converging Nozzle 
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Figure 43. Aluminized Propellant Internal Data; Original 3D Nozzle 


4. Exhaust Data 
Table 14 shows test conditions and results for the 15 successful exhaust meas- 


urements. 






Table 14. ALUMINIZED PROPELLANT EXHAUST RESULTS 
Distance 
from Exit 


Run Num- Nozzle 
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Raw data obtained in runs 9, 10, 11, and 12 were smoothed (maximum of two 
rings) prior to Malvern size distribution calculations to remove effects of an internal re- 
flection caused bv optics on the Malvern protective enclosure. Figure 44 shows the or- 
iginal and corrected scattering data for run number 10. The reflection appeared after 
changing the inline optics and Malvern protective enclosure in an unsuccessful attempt 
to accommodate the shorter cut-off distance of the 63 mm lens. The reflection was 
finally eliminated after approximately ten days of optics adjustments. Data taken during 
this period did not appear to be adverselv affected. once the data were smoothed. 

Exhaust data are shown in Figures 45, 46, 47, and 48. In most cases, the par- 


ticle size distribution was trimodal with peaks at ten and 20 microns as Well as a sub size 
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Figure 44. Internal Reflection Elimination in the Raw Data 


mode. All obscurations were well within the range expected to provide reasonable D,, 
measurements. 

D4, was plotted against chamber pressure, measurement distance from the exit 
plane, and measurement distance from the nozzle throat, but no significant correlation 


was found for D, . 
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rigure 45,  Aluminized Propellant Exhaust Data; Medium Converging Nozzle 
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Aluminized Propellant Exhaust Data; Long Converging Nozzic 
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Figure 45. — Auuninizea Propeliant Exhaust Data; 300 psi Converaing) Diver 
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Figures 49, 50, and 51 show the exhaust D,, plots. The negligible variation in 
D, throughout the plume indicated that the breakup process had been completed at the 
throat. 

All video taped motor firings using the aluminized propellant showed the pres- 
ence of strong shocks behind the nozzle exit. Several shock diamonds were seen behind 
the converging’diverging nozzles, while a single normal shock was visible aft of the 
converging nozzles. For the converging nozzles, data taken immediately aft of the motor 
was inside the shock wave. Particles measured at all other locations were expected to 
show additional breakup effects due to shock wave interaction. None did, further indi- 


cating that the breakup process was completed inside the nozzle. 
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Figure 49. Exhaust SMD as a Function of Chamber Pressure 
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Figure 50. Exhaust SMD as a Function of Distance From Nozzle Exit 
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Figure 51. Exhaust SMD as a Function of Distance From Nozzle Throat 


B. ZIRCONIUM CARBIDE PROPELLANT 
1. Internal 2D Motor Results 
Particle sizing inside the motor with the ZrC propellant was considered more 
accurate due to the lower measured obscurations. (The ZrC propellant contained 1.0 
percent metal; the aluminized propellant metal content was 4.69 percent.) Four internal 
measurements were conducted with the ZrC propellant. Three appear in Figure 52, 
showing particle progression through the nozzle. The fourth was taken inside the long 


converging nozzle, shown in Figure 53. A data summary 1s provided in Table 15. 
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Table 15. INTERNAL ZRC PROPELLANT RESULTS 
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Sauter mean diameters in Figure 52 were corrected using Gulder's empirical 
formulas. The volume distributions presented in the figure explain the decrease in D,, 
through the nozzle. At the nozzle entrance the distribution was trimodal. A slight shift 
towards smaller particles occurred just inside the nozzle, while the distribution data 
taken closest to the throat showed that the larger mode had disappeared. The volume 
distributions for the first two measurements showed no sub size particles, while the data 
nearest the throat bad 2.0 percent of the measured particles below 1.9 microns. 

Unfortunatelv, chamber pressures for the three measurements increased as the 
measurement position approached the throat. Thus, the amount of breakup caused by 


the nozzle could not be separated from the effects of increased chamber pressure. 
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Figure 52. ZrC Propellant Internal Data; Medium Converging Nozzle 
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Figure 53. ZrC Propellant Internal Data; Long Converging Nozzle 
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2. Exhaust Data 
Unlike the aluminized propellant, the ZrC propellant exhibited variation in the 
exhaust gas particle size distribution with nozzle type and beam placement. The nine 


successful exhaust measurements using the ZrC propellant are summarized in Table 16. 






Table 16. ZRC PROPELLANT EXHAUST RESULTS 
Distance 
from Exit 


Run Num- |  Nozzl | | 
B. — ПЕТ! P. (psi) Obscuration D. 
Plane (cm) 
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* measurement position also 1.5 cm off-axis 
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The small amount of data precluded precise determination of the cause of par- 
ticle size variation observed in the plume. The results shown in Figures 54, 55, 56, 57, 
and Figure 58 indicated that particle size in the plume varied somewhat with pressure, 
nozzle contour, and’or exit Mach number. The smallest on-axis D, values were ob- 
tained aft of the converging diverging nozzles, but these measurements Were also subject 
to the highest chamber pressures. 

Youngborg [Ref. 7] used the same ZrC propellant in a small axisymmetric mo- 
tor. In that investigation, it was found that significant breakup occurred within the 
supersonic flow downstream of the throat and that increasing pressure resulted in a 
Siemmereduction in D, at the nozzle throat and at the nozzle exit. These trends are evi- 
dent for the 2D motor, but could not be readily isolated. 

A comparison of distributions for the two motors was also interesting. The 
axisymmetric motor yielded monomodal particles (or monomodal dominated distrib- 
utions) Of approximately 19 microns at the throat. Figures 54 and 55 also show a 


dominant mode at approximately 19 microns. The 2D results showed that a weaker 
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Figure 55. ZrC Propellant Exhaust Data; Long Converging Nozzle 
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Figure 56. ZrC Propellant Exhaust Data; 150 psi Converging/Diverging Nozzle 
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Figure 58. Nozzle Exit SMD as a Function of Chamber Pressure 


mode centered at approximately ten microns also existed. This may have been missed 
in the earlier motor data. At the exhaust of the supersonic flow the axisymmetric data 
showed that some 19 micron particles remained but that three other modes were present: 
a 7-10 micron mode, a 4-5 micron mode, and a mode centered below 2 microns. The 
19 and 7-10 micron modes were particles which did not breakup aft of the throat, 
whereas the smaller particles increased in number, indicating particle breakup in the 
supersonic flow. The data from the two different motors appear to be in reasonably 
good agreement. 

On-axis measurements taken 3.0 cm aft of the motor showed slightly smaller 
D,, values than just aft of the converging nozzle. This could have been caused by shock 
wave interaction. Further experiments are required to isolate the effects of chamber 
pressure, nozzle geometry, and measurement location. 

Two off-axis measurements were made in an initial attempt to determine if any 


size variation occurred across the plume. The two measurements provided ambiguous 
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results: identical beam placement behind different nozzles produced particle sizes that 
in One case Were smaller and in the other case larger than plume centerline values. This 
was not believed to be a nozzle geometry effect. Particle size was expected to be smaller 
at the outer reaches of the plume due to better flow-tracking by small particles. The 
large D,, of 8.1 microns for one of the measurements could have been the result of slag 
shedding from the nozzle, but the volume distribution plot did not show the presence 
of exceptionally large particles. Again, further investigation 1s required to create a sta- 


ustically significant number of samples under the same test conditions. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This investigation produced quantitative particle size information inside the com- 
bustion chamber, inside the exhaust nozzle, and through the exhaust plume under con- 
ditions of varying nozzle geometry and chamber pressure. The two propellants 
examined showed markedly different behavior under similar conditions. 

The aluminized propellant showed no significant variation in plume particle size 
with chamber pressure or nozzle geometry. Nozzle geometry data inside the motor were 
restricted to two different converging sections. Data taken through the 1.05 inch con- 
verging section showed particle breakup as particles approached the throat. The single 
measurement in the 2.00 inch converging section was taken at a chamber pressure sig- 
nificantly higher than that required for comparison with the 1.05 inch converging section 
results. No particle distribution tracking was possible inside the longer nozzle due to 
differences between this nozzle and other designs. Nozzle design had been intended to 
show effects of supersonic expansion on particle breakup. Yet, the aluminized 
propellant exhibited breakup exclusively inside the converging section: size measure- 
ments must therefore be taken inside different nozzles to examine the phenomenon. 
Additional nozzles are required for this investigation. The breakup inside the nozzle 
occurred much sooner than predicted by a rudimentary one-dimensional Weber number 
analvsis. Particle size was small (estimated to be less than eight microns) for all meas- 
urements, including those upstream of the nozzle entrance. This agreed well with com- 
bustion bomb data, which showed little or no agglomeration at 500 psi. 

Particle size distribution was also successfully tracked through the 1.05 inch con- 
verging section using the zirconium carbide propellant. Unlike the aluminized 
propellant, variations in exhaust plume particle size were noted. Sparse measurements 
and unpredictable chamber pressures produced a data set which did not allow separation 
of the effects of chamber pressure and nozzle geometry on particle size. Yet, the ob- 
served effects of pressure and supersonic flow on the dominant modes were in good 
agreement With the axisymmetric motor data reported in Reference 7. 

Larger particles were measured throughout the combustion and expansion process 
using the ZrC propellant than were seen using the aluminized propellant. The larger 
measured size could have resulted from increased agglomeration, slower combustion, or 


reduced particle breakup. 
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Particles produced by either propellant were notably smaller at higher pressures in- 
side the motor. This phenomenon was also observed but not fully verified in exhaust 
measurements with the ZrC propellant. 

Initial attempts at measuring off-axis particle size in the ZrC exhaust plume 
produced ambiguous results, requiring further investigation. 

Throughout the course of this investigation, several factors with potential to in- 
crease test-to-test repeatability were identified. One such factor was propellant slab 
width. Slab thickness varied several hundredths of an inch (between slabs and at differ- 
ent points on the same slab.) This was believed significant enough to cause poor RTV 
sealing or propellant compression inside the motor, resulting in unpredictable chamber 
pressure. Better slab thickness quality assurance is currently not possible at the Naval 
Postgraduate School due to safetv restrictions for cutting propellant. A second im- 
provement could be obtained with. Malvern modification. Laser diffraction particle 
sizers produced bv Malvern Instruments, Inc. are best suited for measurements of un- 
changing phenomena. In the case of the 2600HSD model, this is a software restriction. 
To optinuze effectiveness for the rocket motor environment, a customized software up- 
date from Malvern Instruments could allow many size measurements of short duration 
to be made during each experiment. Current software limits the number of sample 
measurements that can be taken during each experiment due to the excessive time re- 
quired for data storage and processing between measurements. An additional improve- 
ment would be accomplished by extending the range of empirical obscuration 
corrections. Unfortunately, this would entail exhaustive experimentation. 

Improvements in theoretical flow modeling for comparison with experimental results 
are required. The existing breakup theory was not used to its full potential due to lim- 
ited knowledge of high temperature phenomena, such as gas viscosities and droplet sur- 
face tensions. In depth three dimensional flow modeling would also be of benefit. 

Significant particle breakup occurred earlier in the nozzle contraction than predicted 
by one-dimensional calculations. If particle breakup does occur at much lower Weber 
numbers than currently modeled (due to unknown particle and gas properties) then 


higher port velocities may possibly be used to enhance particle breakup. 
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